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Abstract 



Application of superconducting electrical machinery to 
ship propulsion yields light weight and small volume propul- 
sion plants with Increased system efficiency, reliability 
and lower cost. Such propulsion systems might be the only 
feasible answer for high performance ship types. 

A design method for two-pole synchronous machine with 
superconducting field winding Is developed which will survive 
even under the worse fault conditions. Iron shielded machines 
with the effects of higher harmonics and partial winding 
angles are considered. The method Is then adapted to a 
computer program for ease of obtaining minimum volume and 
weight machines. The results are compared with 20,000 HP 
conventional propulsion generator. It Is shown that great 
weight and volume savings are possible. 
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NOMENCLATURE 



Subscripts 



a 

f 
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0 
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d 

q 

sd 

sq 

Symbols 



armature winding 
field winding 
radial 
azimuthal 

inside (for field quantities — at a radius less 
than inner radius) 

inner (when used with a dimension) 

outside (when used with a field quantity) 

outer (when used with a dimension) 

field quantities, generated by a sheet current 

phase a 

phase b 

phase c 

direct axis-armature 
quadrature axis-armature 
direct axis-damper shield 
quadrature axis-damper shield 



B 

B 

B 

d 
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sat 

2 

av 

w 

sh 



magnetic flux density 

saturation flux density of shield material 

mean-squared field experienced by armature 
conductors 

diameter of primary conductor 

dissipation density of shield material (watts/lb) 



open circuit voltage (RMS) 

voltage behind subtransient reactance 

voltage behind subtransient reactance (per unit) 

magnetic field density 

current 

rated armature current (RMS) 
current density 

rated armature-current density 
rated field-current density 
shield structure thickness 
primary shield thickness 
straight-section length 
total machine length 

effective length for calculating armature 
self inductance 

effective length for calculating field-winding 
self -Inductance 

effective length for calculating field-to- 
armature mutual inductance 

effective length for calculating damper self- 
inductance 

effective length for calculating damper-to- 
armature mutual Inductance 

effective length for calculating field-to- 
damper mutual Inductance 

length of armature conductors 

effective length for calculating eddy-current 
losses 

effective length for calculating shield losses 

self inductance, armature phase 
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fd 
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P 



at 



ft 



ed 
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os 
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’^d" 
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y 



self inductance, field 

mutual Inductance, phase-to-phase 

mutual Inductance, armature-to-damper 

mutual Inductance, f ield-to-damper 

magnitude of mutual Inductance, f ield-to-armature 

mutual inductance, fleld-to-phase-a 

magnitude of mutual Inductance, damper-to-phase-a 

number of turns in armature winding 

number of turns in field winding 

number of pole pairs 

rated volt-amperes 

conduction losses 

eddy-current losses 

shield dissipation 

inner radius, shield 

outer radius, shield 

outer radius, damper shield 

rated terminal voltage 

synchronous Impedance (in ohms) 

per-unit synchronous impedance, normalized to 
Internal voltage, 

per-unit synchronous impedance, with as 
base voltage 

per-unit transient reactance (direct axis) 
per-unit subrtanslent reactance 



Rjj^/R 



ao 

fo 
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wf . 



wfe 



wa 



0 



wae 



T 
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4 > 

0) 

w. 



skin depth 

angular displacement 

field-winding Included angle 

P0 ^ 
wf 

armature winding angle 

p0 

wa 

magnitude of maximum radial stress under three 
phase fault from load 

radial stresses 

sheer stresses 

density 

conductivity of armature conductors 
conductivity of damper-shield material 
rotor-axis displacement from axis of phase-a 
angular frequency, also used for deflections 
electrical angular frequency = pw 
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INTRODUCTION 



Advantages of electrical ship propulsion systems were 
realized years ago. Some of them are: arrangement flexib- 

ility, ease of automation and possible Integration to ship 
service power. Major drawbacks were the relatively high 
cost and weight of system components. 

Over the past ten years or so many studies have been 
done on the possible applications of superconductors to 
electrical power equipment, with most effort devoted to 
transformers with decreasing attention given to transmission 
ac machines, and dc machines respectively. The development 
of hlgh-fleld superconducting wire, suitable for fabrication 
Into windings for electric machines, high vacuum seals and 
rotating dewar technology, made applications of supercon- 
ductors possible for rotating electric machines. 

Among the many different types of superconductors , the 
most useful one Is type II, because of Its ability to carry 
a high current density In a high magnetic field, and the 
first major consequence of this upon machine design Is that 
we may break free from the magnetic limitations of Iron. 
Elimination of Iron will reduce the machine weight and 
provide more armature conductors for a given space, thus 
Increasing the power rating per unit volume. If the flux 
density at the armature conductors Is Increased by a factor 
of 3 over that In a conventional machine, the generated emf 
per unit length of conductor Is thus Increased 3 times. If 
the armature diameter and the rotor speed are left unchanged 
Therefore, for the same output power and current, the 
machine length can be reduced by the factor 3> with a 
corresponding reduction In machine volume and weight. 

Application of superconductors to rotating machines Is 
not without problems. First of all, they will operate only 
at extremely low temperatures. Secondly, perhaps the most 
Important one Is that time varying magnetic fields produce 
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losses which are prohibitive at power frequencies. These 
losses are the limiting factor for the application of super- 
conductors to field windings. Therefore, the configuration 
of a superconducting machine must be a superconducting 
stationary or rotating field winding and a normally conduct- 
ing ac armature winding, which may also be rotating or 
stationary. Configuration for rotating armature has Inherent 
problem of transmitting power across sliding contacts. 

With the advent of rotatable liquid helium transfer coupling^, 
a rotating superconducting field winding became practical, 
thus eliminating the necesity for rotating the armature. 
Different designs and configurations have been studied in 

567 

different countries, among them synchronous machines ’ ’ , 

dc machines , homopolar machines ’ , double armature ma- 

1 2 
chines . and the toroidal reciprocating alternator . Some 

of them are already built or under development. 

Application of superconducting electric machinery to 
ship propulsion systems has been studied over approximately 
the past seven years'^ . A superconducting propul- 
sion system is shown in Figure 1 . The required set of 
performance characteristics from a propulsion plant will 
depend mainly on the ship type being considered. Such 
characteristics are tied with the options, control and type 
of electric machinery. For example, bulk carriers which 
run at their full speed 85% of their operating time will 
not require speed control as much as a naval combatant 
vessel, and the option box may be omitted and a simple 
synchronous-synchronous drive in which the propeller speed 
can be controlled by varying prime mover speed, may be 
used. Destroyers, which operate over their full speed range, 
will require more sophisticated controls and the option box 
may be filled with a cyloconverter if a synchronous system 
is used or dc machinery might be more appropriate. It 
should be pointed out that the choice of thruster and prime 
mover will also affect the electric machinery. 
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Superconducting Electric Propulsion System 



Recently some studies have been done on conventional 

12 13 l4 

electric propulsion of bulk carriers ^ ^ . Sixteen 

tankers, container ships and LNG carriers are being built 
or are under design. Conventional electric propulsion 
systems have been selected due to the fact that flexibility 
of propulsion systems may lead to possible savings in the 
volume of the plant. Gas turbi.ne generator sets can be 
placed on the main deck, reducing the inlet and exhaust 
ducting volume and flow pressui-e drops, thus improving gas 
turbine performance. Added to these advantages, if a super- 
conducting drive is used even more savings are possible. 

IS 

It can be shown that the reduction of specific machinery 
weight from 30 to 10 Ibs/SHP will permit an Increase in 
maximum speed from 35 to 50 knots with no sacrifice in 
payload or endurance for a destroyer type ship. Introduction 
of gas turbines has made specific weights somewhat lower, 
but introduced new problem areas. Newly developed ship 
types such as hydrofoils, SES and catamarans which require 
unusual propulsion plant arrangements, suffer greatly with 
conventional propulsion plants. With increased power levels 
required by high performance ship types, superconducting 
electric propulsion plants may be the only feasible answer. 
The purpose of this thesis is to obtain a design method for 
synchronous machines which can be used to realize possible 
advantages . 
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Chapter 1 

DEVELOPMENT OF A MODEL 
.A Model Configuration 

The geometrical model used for this analysis is 
illustrated in Figure 2 and Figure 3> with approximate 
account taken of the effects of end turns as indicated by 

5 

measurements made on the first MIT experimental machine . 

The axial length of one set of armature end turns is taken 
to be equal to the average armature winding diameter in the 
active section divided by the number of pole pairs, and the 
radial thickness of the end turns is assumed equal to twice 
the radial thickness of the armature winding in the active 
region. These relatively large end turns are assumed 
because there is more copper in the active region and more 
volume is needed in which to cross the phase conductors 
over each other. 

The machine has a smooth outer shell made of either 
laminated iron or a highly conducting material, such as 
copper or aluminum, in order to provide a uniform boundary 
condition and confine the magnetic field within the machine. 
In the case of conductive image shield, circulating image 
currents serve to confine the magnetic field. 

There is a thin highly conductive shield between the 

field winding and armature which rotates with the field 

winding. The purpose of this shell is to Intercept 

asynchronous magnetic fields due to load imbalance and 

space harmonics of the armature field and take large fault 

1 X 8 

torques and radial loads under fault conditions ’ . It 

may also intercept thermal radiation from the room tempera- 
ture armature or a separate thermal radiation shield may be 
used. In any case, the essential purpose of this shell is 
to protect the field windings. Since this shell will not 
be strong enough to take large fault loadings, it should 
possibly be supported by a shell made out of a stronger 
material. In this work, the load carrying structure will 
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Figure 2 

Cross Sectional View 
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Side View 




be called the shield structure and the shell v/hich provides 
actual shielding will be called the primary shield. Shown 
in Figures 2 and 3, the primary shield is inside the shield 
structure, but it may very well be outside the shield 
structure. Consequences of these configurations will be 
discussed in a later section. 

The electrical configuration is shown in Figure 4. 
There, are a total of six windings, three phase, one field 
and a direct and a quadrature axis damper winding. The 
electrical parameters of the machine are the coefficients 
of the flux-current relationship which can be expressed in 
matrix form^^’^"^: 
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Figure 4 

Electrical Configuration 
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This may be expressed by: 
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By doing so, it is assumed that all rotor-to-armature 
mutual inductances vary sinusoidally with rotor position. 
This is equivalent to neglecting higher order harmonics. 

It should be emphasized that higher order harmonics in 
phase-to-phase mutual inductances were not neglected in this 
study . 
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IB Field Analysis 

For the purposes of this study, a two dimensional 
representation Is used. The effects of finite length and 
end turn geometry will be Included In an approximate fashion. 
In the straight, active section of the machine, currents 
are parallel to the axis of the machine and effective 
armature and field current densities are related to terminal 
current by 



J 

a 



2 ^a \t 



e R 2(i_x2) 

wae ao 



( 1 ) 






2 If "ft 



Ve 



( 2 ) 



The field and armature winding are assumed to have a 
uniform density between their Inner and outer radii, and 
to fill the wedged shape space within a winding angle. 
Expressions for the magnetic fields generated by these 
currents are derived In references l6 and 17, and are 
summarized In Table 1. Wherever double signs (± or +) 
appear, the upper sign refers to the laminated Iron shield 
case, the lower to the conductive Image shield case. Unless 
otherwise specified, all current densities take space factors 
Into account. 
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•TABLE 1 
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Table 1 

Magnetic Fields (continued) 
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Table 1 

Magnetic Fields (continued) 
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1C Inductance Parameters 



Self and mutual inductances of and between the various 
windings are calculated by integrating flux linkage of the 
fields given in Table 1, over the area of a winding. 

Reference l6 gives such a calculation and the results taken 
from this reference are summarized in Table 2. 

The electrothermal shield has been treated as a 
fourler series of sinusoidally wound damper windings. The 
justification for this is that each space harmonic field 
will excite currents in this shell that are sinusoidally 
distributed in angle. Further more, each space harmonic is 
orthogonal to all but one of the ficticious windings. Each 
winding is characterized by a self inductance and mutual 
Inductances with field and armature windings. The mutual 
inductance between any two of these damper windings is zero. 
Included in the expressions, Nnp is fictitious number of 
turns of damper winding. 

The expressions given in Table 2 have been left with the 
effective length noted differently for each parameter. 

These lengths may not correspond to any physical lengths of 
machine parts, because of the effects of end turns. Further 
more, since end turns of different coils in a machine will 
have different lengths, we may expect the colls to have 
differing effective lengths. Certain assumptions utilized 
when evaluating effective lengths are discussed in section 
D of this chapter. 
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. TABLE 2 
INDUCTANCES 
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ID Machine Rating Equations 

The volt-ampere rating of a machine Is given by 









where are RMS quantities. 

The internal voltage is given by 
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1/2 



( 3 ) 



(4) 



and 



V, 






X sin'll (see Pig. 5) (5) 

cL 




Voltage-current relationship 
30 



where Is synchronous reactance normalized to internal 
voltage 



and 



. "e'V^-ab'Ia 
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They may be combined to give 
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t 



where C is given by 
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and C„ ^ can be obtained from Table 2. 
snp 

Transient reactances are calculated with the assumption 
that : 



a) Field-winding flux is constant 

b) Damper-winding current is zero 



1 6 1 *7 

These reactances are"^ ’ : 

X ’ = X fl— — 



(3) 






where X^ is per unit synchronous reactance 



"d V, 



= X • 



a V. 



(9) 



Subtransient reactances are calculated assuming that 

flux linked by the damper shield is constant. These are 

K 16,17 
given by ’ 






M 2 
s 

(L -L . 

^ a ab 




( 10 ) 



Since the shell is symmetric, direct and quadrature axis 
subtransient reactances are the same. 
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IE Machine Losses 



There are mainly three different losses in synchronous 

V,- 16 

machines : 



a) Current-conduction losses 



P = 



oX 



( 11 ) 



b) Eddy-current losses: 



ec 32 



2 

where is the mean-squared magnetic field seen by the 

armature conductors; 
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where B^^^is given in Table 1. If the first harmonic of 
B^^^ is considered, for a two-pole (one pole pair) iron 
shielded machine: 
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When this is substituted into (13): 
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c) Shield losses: 

If an image shield is used, circulating currents 
will produce losses: 







6 

s 



( 15 ) 



where is B. evaluated at shield inside radius, and 6 

©s © s 

is the skin depth of the shield material. 

If an iron shield is used, eddy-currents and hysterisis 
will produce losses: 

P , = tt(R ^-R ^)0 , p , d , 

sh os s ^■sh sh sh (Id) 

where d^^ is shield dissipation per unit mass of shield 
material. The effective lengths will be discussed in 
Appendix B. 
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IF 



Fault Effects 



Although elimination of iron reduces the machine vol- 
ume and weight, it degrades the transient behaviour of 
superconducting machines. Large fault torques and stresses 
will be Imposed upon the field winding. In addition to 
fault loading, there will be heat dissipation and field 
current rise, which may drive superconducting field winding 
normal. The protection needed will be provided by a suitable 
shield around the field winding. The expressions for fault 
loadings under various fault conditions are derived in 
references 1 and l8. For us the maximum values are import- 
ant in order to design a machine which can withstand such 
loads . 

The maximum radial loads imposed upon the primary 

shield will be under three phase short circuit from load. 

1 1 8 

These are given by ’ : 






0i 



( 17 ) 



where H ^ and H ^ are the fields outside and inside the 
shell respectively after the fault : 
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and H Is the armature field at primary shield radius: 
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and is the field generated by field windings at primary 

shield radius: 
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( 20 ) 



It should be noted that only the first harmonic fields 
have been taken into account. 

When equations (l8), (19) j and (20) are combined with 
(17) and the maximum stress evaluated when wt^='n', 

a,j,= 2yQH^(l+f)[H^(l+f)+H^0]cos^(pe-(t)^) (21) 



The maximum torque will be applied on the primary 

shield under line-to-line short circuit from load. This 

1 1 8 

torque is given by ’ : 
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( 22 ) 



and the maximum value will occur at wtQ=V2,for which 
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In above derivations a factor voltage behind the 

subtransient reactance, was introduced to correct fault 
loads from open circuit to load in an approximate fashion. 
Figure 6 shows a method to calculate e^" from which 



e 



tf 

f 




A 



l+X^" 



siny/)‘ 



(x^- 



cosy/) ' 



( 24 ) 




j (^a~ ^a) 
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Figure 6 

Vector Diagram of Internal, Terminal and 

Voltages 
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Reaction 



Since, primary shield will not be strong enough to 

resist fault loadings, it must be supported by a shield 

structure. Stresses induced under radial and torque loads 

are derived in Appendix A. 

Under fault conditions, the superconducting field 

winding may be driven normal, due to the heat dissipation 

and the field current rise. Perhaps the most important one 

1 8 

is the field current rise in both dc and ac levels. Allow 

able field current rise will depend on the rated operating 

point with respect to superconductor short sample curve. 

Figure 7 illustrates the allowable field current rise for 

a given operating point. Ac component of the field current 

1 8 

rise will be small compared to the dc component and may 

be neglected. Field current will rise exponentially and 

its final value will depend on primary shield time constant, 

1 1 8 

machine reactances and circuit breaker opening time ’ : 

if = i^^[l + ■ (1-e‘^/^d")] (25) 

d 



where i^ and i„ are the Instant and initial values of 
f f o 

field current respectively; and, T, " is related to the 

l8 

primary shield time constant by , 





( 26 ) 



and for a thin primary shield. 



shield time constant 



18 , 16 . 
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(27) 



This time constant will also determine the machine swings 
under transient conditions, but this was not Included in 
this study. 



38 



Flux density 




Figure 7 

Allowable Field Current Rise 
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IG Design Constraints 

A superconducting coll Is characterized by Its short 
sample test curve. At a given operating temperature, the 
relation between flux density and current density which can 
be carried without driving superconductor normal, constitutes 
such a curve. Any chosen operating point should be to the 
left of this curve. In order to provide stability against 
disturbances of electrical, thermal, or mechanical nature 
which may occur, a material of high electrical and thermal 
conductivity Is Included In the conductor. Copper Is 
usually used for this purpose. Hlgh-purlty aluminum has not 
been accepted because of the problem of making a metal- 
lurgical bond with the superconductor although It has higher 

electrical conductivity. Recently a metallurgical bond has 
22 

been made but the brittle nature of It still limits 
aluminum as a stabilizer. Figure 8 shows Nb T1 superconducto 
short sample test data, with the state of the art line given 

In reference 21. To the right of this line small disturb- 

ance may drive superconductor to normal, according to the 
stability criteria given In this reference. Although such 
a line Is quite arbitrary. It defines a region which Is 
experienced with other magnets. It should be pointed out 
that the higher flux densities are possible with Nb^ Sn 
conductors, but they are more expensive than Nb T1 
conductors and cannot be fabricated Into complex windings. 

Another major limitation comes from the saturation 
considerations of the Iron shield used to confine magnetic 

fields Inside the machine volume. The shield may be locally 

saturated or the maximum flux density carried by the shield 
may exceed saturation flux density of the Iron. These 
considerations will be discussed In detail In Appendix B 
when the computer program Is described. 

Armature current density will be determined by the 
losses and cooling method employed. A figure of 3*0 x' lO^' 
amp/m^could be accepted for a nitrogen cooled armature. 
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FIELD CURRENT DENSITY, 10® RMF/M^ 

Figure 8 

Superconductor Short Sample Test Data 
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structural integrity of the machine structure will also 
constitute a problem especially when It Is degraded under 
fault conditions. The most critical part of the machine 
will be the shield structure and must be designed to take 
large fault loadings. This will be discussed in Appendix 
A and B, In detail. 
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Chapter 2 
RESULTS 

2 A Parameters Used In the Design Study 

For the purpose of this study, calculations were 
performed for a 20,000 HP, two-pole synchronous generator. 
Input parameters to the computer program are given In Table 
3. First harmonic field In the bore of the field winding 
Is shown In Figure 9 along with the superconductor short 
sample test curve. Shield structure is assumed to be 
stainless steel and low values of working stresses are 
chosen due to the fact that heat dissipation will degrade 
material properties. The 120 degree field winding angle 
has an advantage of eliminating third harmonic component of 
field thus preventing flow of third harmonic currents in 
delta connected machines. In addition to the input para- 
meters given in Table 3> field winding inside radius of 
2.76 inches and outside radius of 5-0 inches were chosen, 
and the iron shield was assumed to be on the armature end 
turns (t = 0). 

Although the program could be run several times with 
different field winding radii and input parameters to get 
absolute minimum weight and volume machine, this process 
was omitted here and only a design for a given set of input 
parameters defined in Table 3, was obtained. 
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Figure 9 

First Harmonic Field in the Bore of the Field Winding 



TABI.E 3 



DESIGN CONSTANTS 

DESIGN HORSEPOWER= 20^150.25 

MVA RATING= 15.25 

POWER FACTOR=1.00 

DESIGN RPM=3600.0 

RATED FIELD CURRENT DENSITY=1.15 10 8 A/M2 
RATED ARMATURE CURRENT DENSITY=3.00 10 6 A/M2 
MAXIMUM FLUX DENSITY IN FIELD WINDING=5.00 TESLAS 
FIELD WINDING ANGLE=120.00 DEG. 

ARMATURE WINDING ANGLE=60.00 DEG. 

ARMATURE SPACE FACTOR=0.27 

ARMATURE FILAMENT DIAMETER=0 . 04 0 IN. 

YOUNG MODULUS OF SHIELD SUPPORT MATERIAL=29 • 0 10 6 PSI 

POISSON RATIO OF SHIELD SUPPORT MATERIAL=0 . 300 
ELECTRICAL CONDUCTIVITY OF SHIELD=6.0 10 7 MHO/M 
DENSITY OF SHIELD SUPPORT MATERIAL=0 . 290 LBS/IN3 
WORKING STRESS IN TENSION=50000 . 0 PSI. 

WORKING STRESS IN COMPRESSION=50000 . 0 PSI. 

WORKING STRESS IN SHEAR=30000 . 0 PSI. 

CIRCUIT BREAKER OPENNING TIME=10.0 CYCLES 
RATED FIELD CURRENT=680 . 0 AMP. 

ALLOWABLE FIELD CURRENT RISE=420.0 AMP. 



45 



2B Effect of Primary Shield Location 

Figure 10 shows the effect of primary shield location 
where the induced currents flov: under fault conditions, 
for different air gap lengths. Calculations are done for 
the same set of input parameters as defined in Table 3 * The 
case where the primary shield is inside the shield struct- 
ure, produces less radial stresses than the case where 
primary shield is outside the shield structure, under 
three phase short circuit from load, thus may require less 
shield structure and alrgap length and therefore less 
machine volume. Two computer outputs given in Tables 4 
and 5j show this effect. As can be seen, the results are 
not very different from each other. Although the case 
where primary shield is inside shield structure gives 
smaller machine volume and length, the direct exposure of 
the shield structure to the trapped magnetic fields may 
induce higher heat dissipation thus increasing the 
refrigeration power required. 
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Figure 10 

Effect of Primary Shield Location 



• TABLE 4 



FIELD WINDING OUTSIDE RADIUS= 5-000 IN. 

FIELD WINDING INSIDE RADIUS= 2.760 IN. 

ARMATURE INSIDE RADIUS= 6.451 IN. 

ARMATURE OUTSIDE RADIUS= 9-354 IN. 

SHIELD INSIDE RADIUS= 12.257 IN. 

SHIELD OUTSIDE RADIUS= 19-887 IN. 

TOTAL LENGTH=58 .848 IN. 

STRAIGHT SECTION LENGTH= 27-237 IN. 

MACHINE VOLUME= 42.313 FT 3 

WEIGHT OF ACTIVE PARTS= l64l8.903 LBS. 

PRIMARY SHIELD IS OUTSIDE THE SHIELD STRUCTURE 
AIR GAP LENGTH= 1.4514 IN. 

PRIMARY SHIELD OUTSIDE RADIUS= 6.364 IN. 

SHIELD STRUCTURE OUTSIDE RADIUS= 6.364' IN. 

PRIMARY SHIELD THICKNESS= 0.0 IN. 

SHIELD STRUCTURE THICKNESS= 1.3410 IN. 

SHIELD DEFLECTION TOWARDS FIELD WINDING UNDER THREE PHASE 
FAULT= 0.0220 IN. 

SHIELD DEFLECTION TOWARDS ARMATURE WINDING UNDER THREE 
PHASE FAULT=-0.0205 IN. 

SYNCHRONOUSREACTANCE P.U. =0.2533 

TRANSIENT REACTANCE P.U. =0.1875 

SUBTRANSIENT REACTANCE P.U. =0.1283 

MACHINE LOSSES PERCENT OP MACHINE RATING=0.751 

MAX. RADIAL STRESS UNDER THREE PHASE FAULT= 3482.8 PSI 

MAXIMUM TORQUE UNDER LINE TO LINE SHORT CIRCUIT FROM LOAD 

10.303 
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P.U 



■ TABLE 5 



FIELD WINDING OUTSIDE RADIUS= 5-000 IN. 

FIELD WINDING INSIDE RADIUS= 2.760 IN. 

ARMATURE INSIDE RADIUS= 6.306 IN. 

ARMATURE OUTSIDE RADIUS= 9-282 IN. 

SHIELD INSIDE RADIUS= 12.257 IN. 

SHIELD OUTSIDE RADIUS^ 19-887 IN. 

TOTAL LENGTH= 57-950 IN. 

STRAIGHT SECTION LENGTH= 26.773 IN. 

MACHINE VOLUME= 41.667 FT 3 

WEIGHT OF ACTIVE PARTS= I 6169.683 LBS. 

PRIMARY SHIELD IS INSIDE THE SHIELD STRUCTURE 
AIR GAP LENGTH= I .3063 IN. 

PRIMARY SHIELD OUTSIDE RADIUS= 5-026 IN. 

SHIELD STRUCTURE OUTSIDE RADIUS= 6.2941 IN. 

PRIMARY SHIELD THICKNESS= 0.0 IN. 

SHIELD STRUCTURE THICKNESS= 1.2680 IN. 

SHIELD DEFLECTION TOWARDS FIELD WINDING UNDER THREE PHASE 
FAULT= 0.0245 IN. 

SHIELD DEFLECTION TOWARDS ARMATURE WINDING UNDER THREE 
PHASE FAULT=-0 . 0229 IN. 

SYNCHRONOUS REACTANCE P.U. =0.2519 
TRANSIENT REACTANCE P.U. =0.1849 
SUBTRANSIENT REACTANCE P.U. =0.1657 
MACHINE LOSSES PERCENT OP MACHINE RATING=0.749 
MAX. RADIAL STRESS UNDER THREE PHASE PAULT= 3377-9 PSI. 
MAXIMUM TORQUE UNDER LINE TO LINE SHORT CIRCUIT FROM LOAD 
P.U. =8 . 061 
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2C Comparison Bei:ween 20,000 HP Conventional Generator and a 
Superconducting Generator 

The results obtained In previous section are compared 

1 ^ 

with a 20,000 HP conventional propulsion generator . 

Figure 11 shows such a comparison In machine dimensions. 

An additional 15,000 lbs are added to the weight of the 
superconducting generator In order to take account of the 
machine casing, bearings excitation system, etc. Table 6 
lists some of the Information available about the conven- 
tional generator. Unavailability of additional Information 
limits comparison In machine dimensions and weight. As 
far as refrigeration weight and volume are concerned, there 
will not be anything to say unless rotor Is structurally 
and mechanically designed. In addition refrigeration 
system may be Integrated to provide for other superconduct- 
ing electric machinery. 



TABLE 6 

20,000 HP Conventional Generator Characteristics 





Maximum Continuous 


2 Minute 




Rating 


Capacity 


KVA 


15,250 




RPM 


3,600 


900 


Volts 


4,700 


823 


Power Factor 


1.0 




Frequency, Hz 
Amperes 


60 


300 % Max 
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Figure 11 

Comparison Between 20,000 HP Conventional Generator 
and a Superconducting Generator 



Chapter 3 

CONCLUSIONS AND RECOMMENDATIONS 

One major conclusion that can be drawn from the previous 
results is that, superconducting machines can be built to 
withstand fault loadings and still have low weight and 
volume . 

Since most of the weight comes from the iron shield, 
replacement of it by a conductive image shield may lead 
to additional weight savings, although image shield may 
introduce additional losses. 

Replacement of copper by aluminum as a stabilizer will 
also reduce field winding weight if a suitable fabrication 
technique can be found. 

Conversion of the computer program for many pole 
machines will increase its usefulness and provide informa- 
tion about low speed (100-200RPM) synchronous propulsion 
motors . 

Although a simple technique could be employed which 
changes input parameters systematically and gives dihection 
towards minimum weight or minimum volume machine, it is the 
belief of the author that such technique would require a 
great amount of computer time and money. Furthermore, 
minimum weight or volume machine may not be best as far as 
rotor mechanical design and refrigeration input are 
concerned . 

Perhaps, the most questionable assumption is that the 
weight densities for individual machine parts are uniform, 
especially for the rotor structure. Correct estimation of 
rotor weight can be done only by structural and mechanical 
design. Critical speeds of the rotor may not be acceptable 
and in that case rotor structure or length should be 
modified. Rotor lenght should be Increased as much as 
possible in order to increase thermal resistance and thereby 
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reduce heat leak into cold region. Critical speeds will, 
however, limit the rotor length. Inclusion of above consider- 
ations will require longer computation time but may give 
more useful information about superconducting machines. 
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Appendix A 

SHIELD DEFLECTION AND STRESSES UNDER FAULT CONDITIONS 

It has been assumed that only the shield structure will 
be effective for carrying fault loads. Figure A-1 shows 
radial loads Imposed upon simply supported shield structure 
under three phase fault, with the addition of centrifugal 
loading. It is easier to separate these two loadings and 
superimpose them after each solution is obtained separately. 

Reference 19 gives three differential equations for 
deflections : 



32u , 1-v g2u , 1+v 3^u v _ p, 
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2a 3x36 2 3x2 ^ p" 36 12^ 
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+ h- 
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a 3x2 36^ a ^ 3 0 ** 



2-v 3^u 



12 



a 3x236 a 36 






Eh 



(A-lc) 



Simply supported boundary conditions imply that, 

at X = 0 and x = L 
u=0,w=0, N=0,M =0 

X X 

These boundary conditions can be satisfied if one chooses : 



54 




L H 




Figure A-1 

Three Phase Fault Radial and Centrifugal Lo .ding on Shield 

Structure 
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7T X 

u = cos— ^ 2 A cos n 6 
L n 



u = sin^ 2 B sin n0 
L n 



(A-2a) 

(A-2b) 



w = sin-^ 2 C 60 s n 0 
L n 



(A-2c) 



It was shown in Chapter 2-F that the radial stresses 
depend on the angular position along the circumference. 
They are given by equation 21. 



y C) TTX 

cos 0 = 77— (1 + cos 20) x sin— 7 - 
3 ro 2 L(A-3) 



When equations (A-2) and (A-3) are combined with (A- 1 ), 
and the last two terms neglected In equation (A-lb) because 
they will be small for thin shells compared to the other 
t e rms , 



A^[2TT^+(l-v)A^n^]-B^[ (1+v) ATrn]+C^[2vATT] = 0 (A-i|a) 



A^[3 (l+v)TTnA]-B^[ 3 (l-v)7r^+6A^n^ ]+Cj^[6A^n]=0 (A-4b) 

Aj^[vir3A]-B^ nA^ j^3+n^[ (2-v)iT^ + A^n^]j+C^[3A^+r)^ (TT^+A^n^)^] 



where 




12 (l-v^) 



and for n =0 and 2 , Po~^ 2 “^’ result which comes directly 

from equation (A-3) 
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The coefficients obtained by solution of equation (A-4) 
are too complicated to write explicitly. The results are 
summarized in Table A-1. The stresses can be obtained by. 



a = E 



X 






(A-5a) 



Or. = 



[e„+V£^-z(X^+vX^)] 



e l-v2 '-''0 ''^'X ' 0 ' X 



(A-5b) 



’’■X0 2 (1+v) 



(A-5c) 



where 



e = 

X 8x 
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.Q _ 1 ^ w 
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( A-6c ) 



V - W 



( A-6d ) 



Y = ^ ^ ()L + 1 

0 a 30 ^a a 30^ 



( A-6e ) 



Y = 1 rl ^ >1 

x0 a ^2 3x 3x30^^ 



(A-6f ) 



and z is the distance from the mid-surface (positive 
towards the axis of the shell). Equation (A-6) can be 
combined with equation (A-5) to give stresses in axial and 
circumferential directions. 
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TABLE A-l 

Solution of equation A-4 



0 

4 2 

L n cTo 
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6(1+v)ttA-^[7t2+(1-v)2A2] 






(Tr2+4A'')]-2A2[3+nH(2-v>7T2+4A2]]j . [6(1 +v)ttA- 

[7r2+(l-v)2A2]j-j^i^^^li±:^ _3(l_v)Tr2-24A2j . 

[2Tr2+(i-v)4A^] . [3A2+n2(Tr2+4A2)2]/[2vA7rj]j 
B 2 • [2 (1+v) A7t]-A2 . [2tt^ + (1-v )4 A^ ] 



CSA^+n" 

12A 



V7T 



3V7TA- 



2vA7r 



58 



Since the shield will rotate with field winding, 
centrifugal stresses will be Induced with associated 
deflections. The centrifugal loading may be considered as 
loading due to uniform Internal pressure given by: 

pw^.ah 

^ " g~ (A-7) 



Reference 19 gives the solution for a circular cylind- 
rical shell with a uniform internal pressure with ends 
simply supported: 



w 



PL^ 



2 sln°: • slnh°: 
cos2<^ + cosh2“ 



sing (x-^) 



sinhg (x-|^) 



2cos°^ . cosh 
cos2“: + cosh2<^ 



cosg (x-^) 



coshg (x-|!-) 



(A-8) 



where 




and 



<X 



= gL 



2 



Equation (A-8) is the deflection in vertical direction. 
From symmetry, the component v of the displacement In the 
circumferential direction vanishes. The expressions for 
the strain components then become 



e 



X 



du 

dx 



(A-9a) 



w 

= - I 



and it can be shown that 

du _ w 
dx ^ a 



19 



(A-9b) 
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(A-10) 



Equations (A-5), CA-6 ) , ( A-7 ) , (A-8), (A-9), (A-10) 
and Table A-1 can be combined to get total deflections 
and stresses along the shell at any point on circumference. 
The component w of the displacement will be maximum at 
midlength (x=L/2) when 9 = 0 and 9 =tt/ 2 denoted by w^ and ^2 
respectively. These two deflections will determine the 
required airgap along with shield thickness. If the 
material properties are chosen, it is possible to calculate 
stresses at these points, top and bottom surfaces. Since 
the shear stresses will be zero, these are actually prin- 
cipal stresses. Maximum shear stresses can be obtained by. 




Under line-to-line short circuit from load, large fault 

torques are applied on the shield. The magnitude of these 

28 

torques is given by . The shield structure must be strong 

enough to carry these stresses. For the shear stress in the 

23 

shell, according to the maximum shear theory 



(T 



x6 



) 

max 



2 






(A-12) 



in this case ^9 will be assumed to consist of centrifugal 
stresses only, that is. 



2 2 
= £w_ri_ 

9 g 



(A-13) 



and 



x9 



2R T 
so 

F(R ■* R , ' 
^ so - si 



(A-14) 



Equations (A-13), (A-l4), and (A-12) may be combined to give 
maximum shear stresses. 
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Application of this fault torque may cause buckling of 
shield. Reference 20 gives criteria for critical buckling 
stresses . 

• = 0.85 (7rVl2)(l-v2)"5/8E(h/a)5/^(a/L)^/^ 

(A-15a) 

for 10(h/a)^^^<L/a<3(a/h)^^^ 



= 0.272 (1-vM 3 /^E(h/a)^'^^ 

1/p 

for. L/a>3(a/h) ^ 



(A-15b) 
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Appendix B 



COMPUTER PROGRAM LISTING AND DESCRIPTION 



The first harmonic field in the field winding for a 
two pole (p=l) machine is given by (Table 1). 



B, 



2y^ sin 

0 I 



rfw 



wfe 



sin( 9-(J)) 



Stt 



- 2 - 



^fi 



+3 



R 



fo 



-1 



R 



-1 



R 



fo 



R. 



a-y3) 



(A-16) 



The maximum flux density in the r direction can be 
found by taking the derivative of with respect to r 

and setting it equal to zero, which gives a simple result 
for a two pole machine 



r = 



R 



fi 



(A-17) 



and the maximum flux density in the field winding is 



B 



max 



2]Jo sin 


2 


^fo 


37T 



-3y+3+ 


H 


2 3 1 

(1-y^) 




S 





(A-18) 



This linear relation between B and defines the 

max f 

magnet line which was shown in Figs. 8 and 7- If an oper- 
ating point is chosen (given and J^) and field winding 

inside and outside radiuses are specified, equation (A-l8) 
can be solved for R„ , shield inside radius. Although 

O 

fields generated by the armature phases try to oppose the 

fields generated by field v;lnding thus reducing B , there 

in 3.x 

will be additional field concentration due to armature end 
turns and higher space harmonics. Thus, calculating B 

IIlCL ^ 

by equation (A-l8) is a conservative assumption. 
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After the location of the iron shield Is fixed, the 
first thing to be checked Is the saturation limitations. 
Pig. A-2 shows the flux linkages through the Iron shield. 
The total flux carried by the shield, for a two pole 
machine 



-A (’’s) !s «= 



de 



(A-19) 



where the first harmonic component of comes from 

Table 1 and Is given by 



®rfo 



2J^ sin 


®wfe 


COS0 


2 


3tt 






(A-20) 



When equation (A-20) Is evaluated at r = R and substituted 

o 

Into (A-19), Integration gives 



>A-A ‘ '"'o 



wfe 



3-n 






(A-21) 



The limitations on saturation require that 



$ 



A-A 



^ B 



sat 



(■’os-Rs) «s 



(A-22) 



When equations (A-21) and (A-22) are combined, the minimum 

shield outside radius, Is 

’ os 



R 



sin 



wfe 



os 



3tt B 



sat 



R^ (1-y^) + R3 



(A- 23 ) 



where B^^^^ Is the saturation flux density of the Iron 
shield. = 1.5 Weber/m^ was used In the program. 

Local saturation of the Iron shield requires that 






sin 



wfe 



3tt 



■R, 



R 



fo 



R. 



(1-y^) ^ B 



sat (A-24) 
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Figure A-2 

Flux Pattern in Iron Shield 



Algebraic manipulation of equation (A-23) and (A-24) 



gives 



R 



os 



- R, 



(A-25) 



This Is to say that shield thickness must be less than 

or equal to shield inside radius. Since R^ is fixed from 

equation (A-l8), calculated R may not satisfy this 

condition. In this case, program rejects this particular 

design and prints a message saying "NO SOLUTION". Same 

message is also printed when equation (A-l8) gives R less 

s 

than R_ . 
fo 

The maximum alrgap length can be found by assuming 
that armature Inside and outside radiuses coincide, that Is 



A 

max 



s fo 



(A-26) 



where t is the radial distance between shield Inside radius 

and armature end turns. In the program, this A is 

max 

6M 



divided into n steps (n=50 in this program) and in each 
iteration airgap length A is incremented by 
Given A, armature inside radius R . 



R . 



ai 



Rfo * 



(A-27) 



and armature outside radius, R 

’ ao 



R 



ao 




+ R . 



ai 



t)/2.0 



(A-28) 



The program then proceeds to sum up the infinite 
series given in Table A-2 of Chapter 2 for ♦ Then 

synchronous reactance is normalized to Internal voltage 
and x„ as given by equat lon(6 )is calculated. Table A-2 
gives the effective lengths used in this study . Equation 
7 is then utilized to give an expression for the straight 
section length, which results in a second degree equation 
for the straight section length. This equation is too 
complicated to write explicitly and will be omitted here. 
Total machine length is given by 



= £ + 2A£ 



(A-29) 



After all machine dimensions are calculated, the 
program then proceeds to calculate machine reactances. 
Synchronous reactance x^ given by( 9) , transient reactance 
x^* given by(8)can be calculated by summing the series. It 
should be pointed out that up to 15 harmonics were used. 

The user is free to choose number of harmonics to be summed. 

To calculate subtransient reactance x^" , the primary 
shield radius has to be specified. Since for a given air- 
gap length, the range of R^ is known and this can be done 
easily by incrementing R^ in each iteration. For a prev- 
iously calculated airgap length A, R^ is 

(m=50 used in this study) R^ = R^.^ ^ A (A-30) 



m 
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TABLE A-2' 



Effective Lengths 
Armature end turn length 



M = R . + 
ai 



R 



ao 



Effective lengths: armature i _ 

self Inductance ^ 



Mutual inductance, field-to- 
armature 

Mutual Inductance, field-to- 
damper 

Eddy-current loss 

Conduction loss 

Field winding end turn length 



1=1 

m 



Afd - ^ 



I . = Z 
ed 



A = 2AZ 
c 



^^f ^fo ^fi 
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Then x^" can be calculated from(lO). 

Next step Is to calculate required primary shield 
thickness, to give time constant such that, field current 
rise will not be higher than the specified value (SUB- 
ROUTINE TIC). Since the maximum field current rise wiQ.l 
be from 25 



i 



f 

max 




(A-31) 



which depends on . If the maximum rise is less than 

^d' 

the specified rise, then the required shield thickness or 

time constant is set to zero. This is to say that the 

primary shield will be a very thin highly conductive layer 

where induced currents flow. In this study, it was assumed 

that the time constant of shield structure is zero. 

Having specified primary shield thickness, the radial 

stresses o under three-phase fault from load is calculated 
0 

(SUBROUTINE STRES) by(21). 

The most important section of the program is SUBROUTINE 
SHIEL which designs the shield structure according to the 
given working stresses. The designed shield must satisfy 
the following conditions: 



1 ) 



Wi + hh + h + W 2 = A 



where h^ is the primary shield thickness, h is the shield 
structure thickness and w^ and W 2 are the deflections in 
two directions. 



2) ^t “ ^h ” V 2 “ ^ primary shield is 

outside the 
shield struct- 
ure 

R, + h ,p = a -if primary 

shield is inside 
shield structure 
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3 ) 




where are Induced stresses in compression, tension 

and shear respectively and W . are the working 

SIC S 1 U sxs 

stresses in compression, tension, and shear respectively. 



A total of 8 principal stresses are calculated after 
¥^ and ¥2 are found by equation (A-4). These are at 9=0 
and 9='n’/2, on the top and bottom surface. All of the stres- 
ses mentioned above are at the midlength where maximum 
stresses occur. Shield length was assumed to cover straight 
section length plus field winding end turns. After checking 
each of them by the requirements specified above, the 
progrm proceeds to check on requirement 4. If, In the 
course of these calculations, any one of the requirements 
Is not satisfied, computations are repeated by incrementing 
shield structure thickness, until specified alrgap Is filled 
with shield structure. Figure A-3 shows the region to be 
scanned by the program. Since the smallest alrgap length 
will give minimum weight or volume machine, the solution 
will return to the main program If all four of the require- 
ments are satisfied. There may be other solutions at 
bigger alrgap lengths which we will not investigate. 



4) 




stresses Induced under llne- 
to-llne short circuit from 



load 




max 
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Machine volume is calculated by 



V 



TT R 



OS 



z 



t 



and machine weight by 

W=1T(R 2 , . <1. . Ps+’'(Rfo^-Rfl^Ht • Pf 



+ 81 . AMRao^-Rao • Rpi )Pet-^" >*-s 5 • 2 ss 



where 



- armature density 

p - iron shield density 
s 

p^ - rotor density 

p , - end turns density 
0 1 / 



= il500 Kg/m^ 
= 8000 Kg/m^ 
= ^1500 Kg/m^ 
= 4000 Kg/m^ 



p - shield structure denslty=depends on program 



ss 



input 



These average densities are used in the program. 

For selected values of n, m and increment of h, the shield 
structure thickness solution may not converge or may 
converge at a longer air gap length, if these values are 
small. Smaller step sizes will Increase computation time 
and give true convergence, so some sort of a comprlmlse 
should be done. If the solution does not converge then 
program prints a message saying "NO CONVERGENCE". 
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Figure A-3 

Region Scanned by the Program 

Main program and subroutines are listed in Tables A-3j 
A-4, A-5 j and A-6. Table A-6 gives the sample program 
input which is described below. 

Card 1 Format 3P10.2 



Horsepower Rating 


col. 1-10 


Power Factor 


col. 11-20 


Design RPM 


col. 21-30 


Card 2 Format 2F10.2 




Field winding outside radius, 
Inches 


col. 1-10 


Field winding inside radius. 


col. 11-20 



inches 
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Card 3 . Format 3F10.2 



Maximum flux density In field 
winding, Tesla 


col . 


1-10 


Ratedofleld current density 
X 10 , amp/m‘^ 


col . 


11-20 


Field winding angle, 0^^ 
degrees 


col . 


21-30 


Card 4 Format 4F10.2 






Rated armature current density 
X 10“°, amp/m2 


col . 


1-10 


Armature winding angle, 
degrees 


col . 


11-20 


Armature space factor 


col . 


21-30 


Armature filament diameter. 
Inches 


col . 


31-40 


Card 5 Format 7F10.2 






Shield structure material 
young modulus x 10”°, psl 


col . 


1-10 


Shield structure material 
polsson ratio 


col . 


11-20 


Primary shield conductivity 
X lO""^, mho/m 


col . 


21-30 


Shield structure density, Ibs/ln^ 


col . 


31-40 


Working stresses In tension, psl 


col . 


41-50 


Working stresses In compression, 
psl 


col . 


51-60 


Working stresses In shear, psl 


col . 


61-70 


Card 6 Format 3F10.2 






Circuit breaker opening time, 
cycles 


col . 


1-10 


Rated field current, amp 


col . 


11-20 


Allowable field current rise, amp 


col . 


21-30 
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f 







Card 7 



Format II, If 1 is punched in column one, program 
makes calculation for the case where primary 
shield is under shield structure; and, if 2 is 
punched, calculations are performed for the case 
vrhere primary shield is on the shield structure. 

Format F10.2, distance between shield inside 
radius and armature end turns. Inches col. 11-20. 
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CPP=CP*CSNP 
SUM=SUm+CPP 
30 C*'*NTIhUE 

XAP=(AJ«( (RAO/RFO)«»3)«SUM)/(DSQRT(2.0DO)»FI«(DSIN(RTWA/2.0) )*(DSI 
1N(RTWF/2.0) )«(1.0-Y*»3)*C1P) 
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A15=AU«»? 

XOP=XD* (1.0-1.5*A15/(SUM^*SUM1) ) 
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W6=3.U1592653D0* {ARO»«2-ARI»»2)*SL1*39.3700787^DO»GM 
W5=W1+W2 + W3'*’WA 

PCC = 3,0*RTWA*(AJ»*2)«(RAO**2)<^(1,0“X»*2)*SLT/(ASF»6.0007) 
BAV1=1.0/(RAO-RAI) 

BAV2=(2.0»(A.0D-07)»FI*(RFO»»3)»(DSIN(RTWF/2.0) ) » ( 1 .0-Y«*3) /3,0 ) ** 
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904 format { *0» f ‘SHIELD INSIDE RADIL)S=» ♦F7.3*3X* • IN, » ) 



WRITE(6»905) ROSI 

905 FORMAT( »0'» ‘SHIELD OUTSIDE RADIUS= » tF? . 3* 3X* • IN. * ) 
WRITE(6*906) SLTI 
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916 FCKMAT ( *0 •» ‘SYNCHRONOUS REACTANCE P.U.=‘,F6.4) 



917 FOMMAT ( *0» * ‘TRANSIENT REACTANCE P.U.=»fF6.4) 
WRITE(6*9181 XDPP 

918 format ( ‘ 0‘f ‘SUBTRANSIENT REACTANCE P.U,=»*F6.A) 
WRITE(6*919) PPD 
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